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storm  sudden  commencement.  Preliminary  work  on  proton  acceleration  during  the 
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of  the  simulation  code  that  has  been  developed  is  the  ability  to  post-process  data 
runs  after  the  trajectories  of  several  hundred  thousand  ions  or  electrons  have  been 
followed  through  the  inbound  and  outbound  propagation  time  of  the  wave  fields, 
which  is  about  100  seconds  of  real  time.  As  a  first  step  in  the  development  of  the 
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(a)  Objectives 

The  objective  of  this  research  in  magnetospheric  physics  is 
to  advance  understanding  of  the  interaction  of  ring  current  ions 
with  magnetospheric  hydromagnetic  waves.  In  particular, 
excitation  of  westward-propagating  hydromagnetic  waves  by  a  drift- 
bounce  resonance  with  100-200  keV  ring  current  protons  is  being 
investigated  along  with  the  self-consistent  nonlinear  evolution  of 
the  waves  and  the  ring  current  ion  populations  (H+  and  0+)  .  An 
extension  of  the  ring  current  particle  trajectory  tracing  code  has 
been  used  to  study  the  rapid  formation  of  new  electron  and  proton 
radiation  belts,  as  observed  on  March  24,  1991. 


The  methods  employed  in  this  research  program  include  (i) 
numerical  solutions  of  linear  gyrokinetic  eigenmode  equations 
using  existing  computer  codes,  (ii)  simulations  of  the  motion  of 
ring  current  ions  in  hydromagnetic  waves  using  an  existing 
gyroaveraged  Hamiltonian  particle  code,  (iii)  development  of  a  3D 
dipole  geometry  magnetohydrodynamic  (MHD)  code  and  (iv)  coupling 
of  the  test-particle  and  MHD  codes  to  produce  a  hybrid  MHD- 
gyrokinetic  code.  The  results  will  provide  the  first  simulation 
of  the  excitation  and  nonlinear  evolution  of  magnetospheric 
hydromagnetic  waves  due  to  ring  current  ions,  and  of  the  self- 
consistent  effect  of  the  waves  on  the  ions.  The  reduction  of 
storm  time  ring  current  fluxes  back  tc  pre-storm  levels  affects 
the  background  magnetic  field  and  the  higher-energy  radiation  belt 
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spacecraft  data  will  be  made  to  test  the  theoretical  predictions 
and  yield  new  theoretical  insights  into  observed  wave  properties 
and  associated  particle  measurements. 

(b)  Accomplishments 

Radiation  Belt  Results 

Major  progress  was  made  during  this  first  year  in  obtaining 
an  explanation  for  the  rapid  formation  of  new  electron  and  proton 
radiation  belts  during  the  March  24,  1991  SSC.  A  relativistic 
generalization  was  implemented  of  the  guiding  center  test  particle 
code  developed  by  Li  et  al.  (1993a)  to  follow  the  trajectories  of 
ring  current  ions  interacting  with  hydromagnetic  waves  in  the 
inner  magnetosphere.  The  new  code  was  used  with  a  model  for  the 
incoming  and  outgoing  magnetosonic  wave  imposed  on  the 
magnetosphere  by  the  SSC  shock  impact  on  March  24,  1991.  The 
acceleration  of  pre-existing  trapped  electrons  from  L  =  4  -  9  by 
the  incoming  pulse,  due  to  the  induction  electric  field  associated 
with  the  magnetospheric  compression,  has  been  described  in  the 
attached  reprint  (Li  et  al.,  1993b).  Remarkable  agreement  was 
obtained  with  the  electron  drift  echoes  reported  by  Blake  et  al. 
(1992)  at  the  L  =  2.5  nightside  position  of  the  CRRES  satellite 
following  the  SSC,  see  Fig.  1  of  the  attached  reprint. 

Preliminary  work  on  proton  acceleration  during  the  same  event 
was  reported  at  Fall  AGU  (Kotelnikov  et  al.,  1993).  The  proton 
source  population  was  taken  to  be  the  solar  wind  shock  accelerated 
protons  whose  flux  was  observed  to  increase  in  the  Aerospace  and 
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other  instruments  on  CRRES  preceedi.ig  the  SSC  (Blake,  Gussenhoven, 
private  communication;  Shea  and  Smart,  1993) .  A  softer  power 
spectral  index  e~N  was  assumed  for  the  protons  (N  =  2)  than  for  the 

electrons  (N  =  8)  ,  consistent  with  CRRES  measurements  and  those  by 
a  Japanese  satellite  that  was  on  the  dayside  at  the  time  that  the 
SSC  struck  the  magnetosphere  (Blake,  private  communication) .  The 
same  model  wave  field  proved  to  be  very  efficient  at  producing 
proton  drift  echoes  in  our  simulations  with  drift  periods  in 
agreement  with  Aerospace  omnidirectional  data. 

An  important  feature  of  our  simulation  code  is  the  ability  to 
post-process  runs  after  we  have  followed  the  trajectories  of 
several  hundred  thousand  ions  or  electrons  through  the  inbound  and 
outbound  propagation  time  of  the  wave  fields  (~  100  s  real  time) . 
In  the  post-processing  we  can  vary  the  assumed  initial  spectral 
index  choosing  a  subset  of  the  particle  trajectories  followed,  and 
incorporate  detector  response  and  spacecraft  position  in  order  to 
obtain  a  best  fit  of  our  simulated  particle  fluxes  to  those 
measured  by  CRRES  instruments.  This  provides  us  with  information 
about  the  source  population  without  repeating  runs  and  will  allow 
us  to  make  further  comparisons  with  measurements  made  by  different 
instruments  on  CRRES.  We  particulary  look  forward  to  examining 
the  PROTEL  data  for  double  peaking  in  the  ion  drift  echoes  evident 
in  the  omnidirectional  data  (Blake  et  al.,  1992)  .  An  explanation 
for  this  double  peaking  was  suggested  in  the  attached  reprint,  and 
our  ability  to  reproduce  such  features  for  an  instrument  with 
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narrower  energy  and  pitch  angle  resolution  will  place  additional 
constraints  on  our  assumed  source  population  and  field  model. 

Rina  Current  Results 

As  a  first  step  in  the  development  of  the  three-dimensional 
MHD  fluid-gyrokinetic  particle  hybrid  code,  we  have  developed  a 
three-dimensional  MHD  version  of  the  code  in  generalized 
orthogonal  coordinates.  We  are  currently  testing  this  code.  The 
code  has  been  tested  with  one,  two,  and  three-dimensional 
Cartesian  coordinates  by  successfully  reproducing  the  MHD  waves 
and  verifying  energy  conservation.  We  are  now  testing  the  code 
with  curvilinear  coordinates,  specifically  with  dipole 
coordinates.  First  we  tested  a  two-dimensional  version  of  the 
code  in  dipole  coordinates  successfully  with  an  oscillating 
magnetic  flux  bundle  isolated  from  the  boundaries  of  the 
simulation  domain.  The  results  obtained  using  the  dipole 
coordinates  were  equivalent  to  those  using  the  Cartesian 
coordinate  code.  Then  boundary  conditions  appropriate  for  the 
magnetosphere  were  introduced.  Specifically,  a  perfectly 
conducting  boundary  is  implemented  for  the  ionospheric  boundary,  a 
constant  value  or  constant  gradient  boundary  is  implemented  for 
the  radial  boundary,  and  a  periodic  boundary  is  assumed  for  the 
azimuthal  direction.  We  have  reproduced  a  toroidal  field  line 
oscillation  at  roughly  the  WKB  field  line  resonance  period  using 
the  two-dimensional  version  of  the  code.  Further  work  is 
continuing  on  the  boundary  conditions  and  testing  of  the  code  in 
three-dimensional  dipole  geometry. 
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Several  new  results  have  been  obtained  from  numerical 
solutions  of  the  linear  gyrokinetic  eigenmode  equations.  These 
results  are  described  in  the  paper  "Anisotropic  Alfven-Ballooning 
Modes  in  the  Earth's  Magnetosphere"  by  Chan,  Xia  and  Chen, 
accepted  in  November  1993  for  publication  in  the  Journal  of 
Geophysical  Research.  In  this  paper,  motivated  by  spacecraft 
observations  which  show  that  the  perpendicular  plasma  pressure  Px 
is  typically  larger  than  the  parallel  pressure  P||,  we  focused  on 

the  previously-neglected  effects  of  finite  pressure  anisotropy . 
Neglecting  all  kinetic  effects  and  using  a  self-consistent  MHD 
equilibrium  obtained  from  a  multiscale  perturbation  expansion  of 
an  anisotropic  Grad-Shafranov  equation,  we  numerically  solved  the 
corresponding  eigenmode  equations  using  shooting  methods.  The 
main  results  are:  (1)  The  field  line  eigenfrequency  can  be 
significantly  lowered  by  finite  pressure  effects.  (2)  The 
parallel  mode  structure  of  the  compressional  magnetic  component 
can  become  highly  peaked  near  the  magnetic  equator  due  to  large 
perpendicular  pressure.  (3)  For  the  isotropic  case  P||=  Px  =  P 

ballooning  instability  can  occur  when  the  ratio  of  the  plasma 
pressure  to  the  magnetic  pressure,  |3  =  P/(B2/8tc),  exceeds  a 

critical  value  of  about  3.5  at  the  equator.  (4)  Compared  to  the 
isotropic  case  the  critical  beta  value  is  lowered  by  anisotropy. 

(5)  We  use  a  "(3-8  stability  diagram"  to  display  the  regions  of 

instability  with  respect  to  the  equatorial  values  of  the 
parameters  p  and  5,  where  p  =  (1/3)  ({3||  +  2|3X)  is  an  average  beta 

value  and  8=1-  Pjj/Px  is  a  measure  of  the  plasma  anisotropy. 
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The  diagram  is  divided  into  regions  corresponding  to  the 
firehose,  mirror  and  ballooning  instabilities.  (6)  It  appears 
that  observed  values  of  the  plasma  pressure  are  below  the  critical 
value  for  the  isotropic  ballooning  instability  but  it  may  be 
possible  to  approach  a  ballooning-mirror  instability  when 

Pi/P||  >  2. 
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Progress  on  ring  current  and  radiation  belt  code 
development 

X.  Li  (Same  title  as  GRL) 

UC  Berkeley 
Dartmouth 

National  Research  Council  of  Canada,  Herzberg  Institute, 
Ottawa 

Philips  Laboratory,  Bedford,  MA 

(ii)  Consultative  and  advisory  functions 

Denton  visited  Phillips  laboratory  in  June  and  presented  a 
progress  report  on  ring  current  code  development  and 
radiation  belt  studies.  Li  visited  Philips  Laboratory  in 
December  and  consulted  with  Mullen  and  Gussenhoven  on  CRRES 
low  energy  electron  and  proton  data.  Hudson  plans  a  visit  in 
January  to  discuss  application  of  our  model  to  CRRES  PROTEL 
data.  Hudson  has  been  in  frequent  contact  with  Blake  of  The 
Aerospace  Corporation  on  more  energetic  particle  data  from 
CRRES,  and  arranged  a  working  meeting  at  the  Fall  AGU  of  all 
parties  involved  with  the  electron  and  proton  drift  echo 
studies  (GRL  co-authors)  . 
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SIMULATION  OF  THE  PROMPT  ENERGIZATION  AND  TRANSPORT  OF  RADIATION  BELT 
PARTICLES  DURING  THE  MARCH  24,  1991  SSC 

Xinlin  Li1,  I.  Roth2,  M.  Temerin2,  J.  R.  Wygant2,  M.  K.  Hudson1,  J.  B.  Blake3 


Abstract.  We  model  the  rapid  (~  1  min)  formation  of 
a  new  electron  radiation  belt  at  L  ~  2.5  that  resulted 
from  the  Storm  Sudden  Commencement  (SSC)  of  March 
24,  1991  as  observed  by  the  CRRES  satellite.  Guided  by 
the  observed  electric  and  magnetic  fields,  we  represent  the 
time-dependent  magnetospheric  electric  field  during  the 
SSC  by  an  asymmetric  bipolar  pulse  that  is  associated  with 
the  compression  and  relaxation  of  the  Earth’s  magnetic 
field.  We  follow  the  electrons  using  a  relativistic  guiding 
center  code.  The  test-particle  simulations  show  that  elec¬ 
trons  with  energies  of  a  few  MeV  at  L  >  6  were  energized 
up  to  40  MeV  and  transported  to  L  ~  2.5  during  a  fraction 
of  their  drift  period.  The  energization  process  conserves 
the  first  adiabatic  invariant  and  is  enhanced  due  to  reso¬ 
nance  of  the  electron  drift  motion  with  the  time- varying 
electric  field.  Our  simulation  results,  with  an  initial  W~* 
energy  flux  spectra,  reproduce  the  observed  electron  drift 
echoes  and  show  that  the  interplanetary  shock  impacted 
the  magnetosphere  between  1500  and  1800  MLT. 

Introduction 

At  3:41  UT  on  March  24,  1991  an  almost  instantaneous 
formation  of  a  new  electron  radiation  belt  was  observed  by 
the  CRRES  satellite  which  was  fortuitously  situated  at  a 
distance  of  2.55  Re  near  the  equatorial  plane  at  a  local 
time  of  0300  MLT  [Vampola  and  Korth,  1992;  Blake  et 
al.,  1992].  The  creation  of  the  radiation  belt  was  accom¬ 
panied  by  an  electron  and  ion  drift-echo  event  from  which 
Blake  et  al.  [1992]  deduced  that  the  injected  electrons  were 
at  energies  predominately  above  15  MeV  with  an  energy 
spectral  index  of  -6.  This  radiation  belt  lasted  beyond  the 
end  of  the  CRRES  mission  six  months  later.  The  event 
was  apparently  initiated  by  an  interplanetary  shock  that 
compressed  the  magnetosphere  inside  geostationary  orbit. 
The  shock  was  also  detected  by  Ulysses  at  a  distance  of 
2.5  AU  and  by  interplanetary  scintillation  measurements 
from  which  it  was  inferred  that  the  shock  approached  the 
Earth  in  the  afternoon  sector  [Woan  and  Hewish,  private 
communication] . 

The  measured  electric  field  was  seen  as  a  bipolar  pulse, 
predominantly  in  the  azimuthal  direction,  with  a  peak-to- 
peak  magnitude  of  80  mV/m  and  in  the  magnetic  field 
as  an  140  nT  monopolar  pulse  with  a  duration  of  120  s 
[Wygant  et  al.,  1993],  consistent  with  ground  magnetome¬ 
ter  measurements  of  a  narrow  impulse  [Yumoto  and  Sh- 
iokawa,  private  communication].  Wygant  et  al.  [1993] 
argue,  based  on  previous  measurements  of  the  local  time 
dependence  of  field  perturbations  associated  with  SSCs 
[Schmidt  and  Pederson,  1987],  that  the  field  enhancements 
were  much  larger  on  the  dayside. 

The  observed  creation  of  this  electron  belt  is  an  unique 
event  in  space  physics  which  raises  new  questions  about 
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the  formation  of  radiation  belts  in  general  and  about  the 
specific  mechanisms  that  produced  this  electron  belt  in  par¬ 
ticular  [Blake  et  al.,  1992].  In  this  Letter  we  present  results 
from  a  test-particle  simulation  of  the  interaction  of  ener¬ 
getic  electrons  with  electric  and  magnetic  fields  during  the 
SSC.  The  left  column  of  Figure  1  shows  the  electron  data 
from  Blake  et  al.  [1992]  together  with  the  electric  and 
magnetic  fields  from  Wygant  et  al.  [1993].  Questions  that 
these  data  present  and  that  we  address  are:  (a)  How  did 
the  >  6  MeV  electron  flux  rise  in  10  s  at  L  =  2.5  by  four 
orders  of  magnitude?  (b)  What  is  the  source  of  the  elec¬ 
trons  seen  at  L  ~  2.5?  (c)  What  causes  the  oscillations  in 
the  subsequent  electron  flux? 

Model 


We  model  the  unperturbed  magnetosphere  as  a  dipole 
magnetic  field  Be.  The  interaction  of  the  interplanetary 
shock  with  the  geomagnetic  field  compresses  the  magne¬ 
tosphere.  The  compression  is  modelled  as  a  time  depen¬ 
dent  Gaussian  pulse  with  a  purely  azimuthal  electric  field 
component  that  propagates  radially  inward  at  a  constant 
velocity,  decreases  away  from  the  impact  point,  and  is  par¬ 
tially  reflected  near  the  surface  of  the  Earth.  The  modelled 
shock  encounters  the  magnetosphere  first  at  longitude  <j> 0 
and  subsequently  at  other  longitudes.  Explicitly,  in  the 
usual  spherical  coordinates  (r,  0,<b),  where  r  >  Re  is  mea¬ 
sured  from  the  center  of  the  Earth,  <t>  =  0°  at  noon  local 
time,  positive  eastward,  0  —  0°  at  the  north  pole,  the  elec¬ 
tric  field  is  given  by 


E, 


.  =  -e*£0(l+ci  cos(<t>-<j>o)J  jexp(-£2)-c2exp(-v2) 


.  (1) 


where  the  terms  in  the  square  brackets  are  associated  with 
the  compression  and  relaxation  of  the  magnetosphere.  In 
(1)  f  =  [r  +  «o(<  -  trh)]/d.  r/  =  [r  -  v0(t  -  tpk  +  td)}/d,  v0 
is  the  pulse  propagation  speed,  and  d  is  the  width  of  the 
pulse;  C] (>  0)  describes  the  local  time  dependence  of  the 
electric  field  amplitude,  which  is  largest  at  <t>0\  tp/,  =  t,  + 
(c3Re/uo)[1  —  cos(<A-<fo)]  represents  the  delay  of  the  pulse 
from  <f> o  to  other  local  times;  c3  determines  the  magnitude 
of  the  delay;  c2  determines  the  partial  reflection  of  the 
pulse;  Id  —  2.06Re/vq  indicates  that  the  reflection  occurs 
at  r  =  1.03-Re;  and  t,  determines  the  location  of  the  pulse 
at  the  start  of  the  simulation.  In  this  Letter  we  present 
results  with  Eo=  240  mV/m,  cj  =  0.8,  c2  =  0.8,  c3  =  8.0, 
t>o=  2000  km/s,  (,=  80  s,  =  45°  and  d=  30,000  km.  At 
t  =  0  the  pulse  is  at  about  25 Re- 

The  perturbed  magnetic  field  B,„  is  obtained  from  Fara¬ 
day’s  law  and  satisfies  v  ■  (Be  +  Bu,)=0  and  E„,  •  (Be  + 
Bw)=0.  The  model  thus  describes  the  propagation  of  a 
magnetosonic  pulse  through  the  magnetosphere  ignoring 
the  variable  pulse  velocity  due  to  changes  in  density,  tem¬ 
perature  and  magnetic  field.  The  bottom  right  panels  of 
Figure  1  show  the  electric  and  magnetic  field  given  by  the 
model  at  L  =  2.5  and  at  <j>  —  4>0  =  180°  (0300  MLT).  At 
other  local  times  the  form  of  the  pulse  is  the  same  but 
its  amplitude  is  larger.  At  larger  L  values  the  pulse  is 
broader  and  somewhat  larger  since  the  incoming  and  re¬ 
flected  pulses  are  further  separated  and  there  is  less  de¬ 
structive  interference  between  them.  The  model  parame¬ 
ters  were  chosen  to  approximate  the  data  on  the  basis  of 
physical  reasoning  and  some  limited  experimentation.  We 
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Fig.  1 .  a)  Data  from  the  CF.RES  satellite  at  the  time  of 
the  March  24,  1991  SSC  .op  panel  shows  count  rates  as 
a  function  of  time  from  four  energetic  electron  channels 
measuring  integral  counts  above  a  threshold  energy  indi¬ 
cated,  and  also  over  10-50  MeV  [Blake et  al.,  1992].  Middle 
and  bottom  panels  show  the  measured  electric  field  Ey  in  a 
co-rotational  frame  and  the  Bz  magnetic  field  component 
with  a  model  magnetic  field  subtracted,  in  GSE  coordi¬ 
nates  over  the  same  time  interval  [Wygant  et  al.,  1993].  b) 
Simulated  results  in  the  same  format  as  (a)  measured  at  a 
spatial  location  corresponding  to  the  trajectory  of  the  CR¬ 
RES  satellite,  including  the  Aerospace  detector  responses 
in  the  top  panel.  Time  is  measured  from  t  =0  as  deter¬ 
mined  by  (1). 


varied  the  velocity,  v0,  between  750  and  2000  km/s,  the 
field  amplitude,  E0,  between  120  and  240  mV/m,  and  used 
impact  angles  of  0°,  45°,  and  90°.  We  found  that  the  chosen 
parameters  approximate  the  electron  data  well. 

We  follow  the  particles  using  a  relativistic  guiding  center 
approximation  with  uy  =  0,  Er  =  Eg  =  0  [Northrop,  1963] 

W  =  qR1-Ew  +  (2) 

7  eft 

Rx  =  J  x  (-cEw  +  (3) 

where  Rx  describes  the  guiding  center  motion  perpendic¬ 
ular  to  the  instantaneous  magnetic  field  B  =  Be  +  B„,, 
ei  =  B /B  is  a  unit  vector  along  B,  7  =  (W  +  m0?)/moC2 
is  the  relativistic  energy  factor,  W  is  the  particle’s  kinetic 
energy,  Mr  =  pl/2m0B  is  the  relativistic  adiabatic  invari¬ 
ant  and  px  is  tne  particle’s  perpendicular  momentum. 

In  the  simulation  we  follow  336,720  electrons  as  they 
drift  in  the  combined  fields  of  the  Earth’s  dipole  and  the 
modelled  pulse  field,  recording  their  energy,  arrival  time, 
radial  distance,  and  initial  conditions  as  they  pass  vari¬ 
ous  local  times.  We  model  the  relativistic  population  of 
electrons  before  the  event  by  distributing  test-particle  elec¬ 
trons  in  the  equatorial  plane  at  distances  from  L= 3  to  L=9 
in  increments  of  0.1  Re,  in  azimuth  every  3°,  and  at  en¬ 
ergies  between  1  and  9  MeV  in  increments  of  5%.  We  re¬ 
strict  ourselves  to  90°  pitch  angles  following  the  indication 
from  the  particle  detectors  of  an  equatorially  trapped  elec¬ 
tron  population  [Blake  et  al.,  1992].  In  the  post-processing 
stage  each  electron  is  given  a  weight  which  depends  on  the 
assumed  initial  distribution  in  energy  (power  law)  and  in 
L.  The  initial  distribution  in  L  is  modelled  by  a  parabola 
G(£)  =  l—(L—LQ)2/al.  Here,  we  have  taken  the  power  law 


index  -8,  with  the  variation  in  L  determined  by  L0  =  10 
and  a0  =  7.5.  Thus,  given  an  initial  electron  distribu¬ 
tion,  we  obtain  electron  fluxes  and  distributions  at  any 
location  and  compare  the  simulation  results  with  the  C’R- 
RES  observations,  incorporating  the  energy  response  of  the 
Aerospace  energetic  electron  detectors  [Blake  et  al.,  1992], 

Results  and  Discussion 

The  right  column  of  Figure  1  shows  the  simulated  elec¬ 
tron  flux  and  fields  near  the  CRRES  position  at  0300  MLT 
and  L  =  2.5.  Each  point  in  the  upper  panel  shows  the  flux 
over  a  three  second  interval  with  the  corresponding  detec¬ 
tor  characteristics  folded  in  and  integrated  over  L  values  ± 
0.1  Re  around  the  CRRES  position.  The  simulated  L  val¬ 
ues  were  chosen  to  match  those  of  the  satellite  at  the  time 
of  the  first  pulse  and  to  follow  the  orbit  of  CRRES.  The 
agreement  between  the  simulation  and  the  data  is  quite  re¬ 
markable.  This  agreement  includes  the  relative  amplitudes 
and  periods  of  the  first  three  drift  echoes.  Repeated  pulses 
reflect  the  drift  period,  the  sawtooth  shape  reflects  the  fact 
that  the  more  energetic  electrons  which  have  a  lower  flux 
reach  the  detector  sooner,  and  the  abrupt  cutoff  is  due  to 
the  low  flux  of  lower  energy  electrons. 

Since  higher  energy  electrons  drift  faster,  the  magneto¬ 
sphere  acts  as  a  giant  energy  spectrometer.  The  temporal 
slope  of  the  flux  at  a  specific  location  depends  on  the  ve¬ 
locity  dispersion.  From  the  slope  of  the  second  drift-echo 
peak,  Blake  et  al.  [1992]  were  able  to  deduce  that  the  en¬ 
ergy  spectrum  of  the  injected  electrons  on  CRRES  peaked 
at  about  15  MeV,  extended  up  to  50  MeV,  and  had  a  power 
law  spectrum  of  W~6.  In  the  simulations  we  found  that  an 
initial  spectrum  of  W~8  between  1  MeV  and  9  MeV  pro¬ 
duced  the  best  agreement  with  the  data. 

The  magnetospheric  source  of  the  drift  echo  electrons 
can  be  determined  by  comparing  the  measured  electron 
spectrum  in  the  outer  magnetosphere  with  the  measured 
flux  at  L— 2.5,  assuming  that  the  phase  space  density  of 
electrons  is  conserved.  This  implies  that  the  flux  varies 
as  j  oc  L~3  [Lyons  and  Williams,  1984].  Thus  the  mea¬ 
sured  flux  at  L= 2.5  at  14  MeV  requires  an  initial  flux  four 
times  less  at  L= 4  at  6.9  MeV  or  27  times  less  at  L= 7.5 
at  2.7  MeV.  (We  note  that  conservation  of  the  first  adi¬ 
abatic  invariant  for  90°  pitch  angle  relativistic  electrons 
implies  that  the  electron  energy  W  oc  B1^2  oc  L~3?2  since 
W  ~  pc.)  CRRES  observed  the  background  electron  flux 
prior  to  the  drift-echo  event  during  the  inbound  part  of  its 
orbit  from  L= 6.8.  The  CRRES  Cerenkov  counters  provide 
only  an  upper  limit  due  to  contamination  from  penetrat¬ 
ing  solar  protons  that  preceded  the  shock  by  about  a  day 
and  populated  the  magnetosphere  beyond  L  >  4.  This  up¬ 
per  limit  shows  that  there  was  insufficient  electron  flux  at 
L  <  5.  Between  L=b  and  L= 6.8  the  Lockheed  spectrom¬ 
eter  measured  electrons  up  to  5  MeV,  but  was  subject  to 
the  same  contamination  above  2  MeV.  Extrapolating  the 
measured  flux  below  2  MeV  [Nightingale,  private  commu¬ 
nication]  shows  that  there  was  sufficient  flux  beyond  L= 6. 
However,  the  extrapolation  should  be  interpreted  as  an  up¬ 
per  limit  since  it  uses  the  spectral  index  -2.7  to  -4  measured 
below  2  MeV,  rather  than  the  spectral  index  -6  to  -8  in¬ 
ferred  from  the  data  and  the  simulation  at  higher  energies. 
Thus  we  conclude  that  the  source  of  the  drift-echo  elec¬ 
trons  was  beyond  L= 6,  consistent  with  our  simulation  in 
which  most  of  the  flux  is  from  7 <  L  <9. 

The  local  time  of  the  shock  impact  determines  the  width 
of  the  first  drift-echo  peak.  In  Figure  1  the  shock  impacted 
at  1500  MLT.  When  we  choose  1200  MLT,  the  first  peak  in 
the  electron  flux  was  noticeably  wider  than  in  the  data  due 
to  increased  energy  dispersion  of  the  electrons  as  they  drift 
east  from  a  larger  distance.  When  we  choose  1800  MLT, 
the  first  peak  was  noticeably  narrower  due  to  decreased 
dispersion.  For  1500  MLT,  as  shown  in  Figure  1,  the  first 
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peak  is  just  slightly  wider  than  in  the  data.  Thus  our 
simulations  are  in  agreement  with  Blake  et  al.  [1992j,  who 
deduced  from  the  ion  and  electron  data  that  the  shock 
impacted  in  the  afternoon  sector  at  about  1600  MLT. 

Figure  2  shows  the  electron  energy  spectra  at  several 
L  values  integrated  over  a  drift  period.  A  remarkable  as¬ 
pect  of  the  spectra  is  the  sharp  rise  in  the  flux  at  higher 
energies.  As  the  shock-induced  distortion  of  the  magneto¬ 
sphere  propagates  away  from  the  point  of  initial  impact, 
higher  energy  electrons,  which  gradient-drift  faster,  stay 
in  phase  longer  with  the  perturbation  and  thus  are  con- 
vected  in  further.  However,  as  seen  in  our  simulation,  the 
most  energetic  electrons  can  also  outrun  the  optimum  drift 
speed  and  thus  are  not  convected  as  far  inward.  From  Fig¬ 
ure  2  one  sees  that  the  low  energy  cutoff  of  the  electron 
flux  decreases  with  increasing  L.  At  the  lowest  energies 
there  is  again  a  rise  in  the  electron  flux  because  electrons 
originally  at  low  L  values  (3  <  L  <  4)  E  x  B  drift  inward 
without  optimally  matching  the  drift-resonance  condition. 
These  electrons  have,  however,  a  small  phase-space  density 
at  higher  energies  because  they  have  gained  relatively  little 
energy  in  moving  inwards  and  because  of  the  initial  W~8 
energy  spectrum.  Conversely  the  largest  flux  is  contributed 
by  electrons  from  the  largest  L  values  in  the  simulation. 
Most  of  the  simulated  flux  at  L  =  2.5  is  due  to  electrons 
originally  at  L  values  between  7  and  9. 

Low  energy  electrons  are  affected  mainly  by  the  ExB 
drift  and  their  motion  is  indicative  of  the  relative  motion  of 
the  magnetospheric  plasma  and  of  the  magnetopause.  At 
1500  MLT,  the  impact  location  of  the  shock,  zero  energy 
electrons  from  L  =8.5  were  convected  inward  to  L  =  3.84 
and  then  back  to  L  =  5.51  by  the  reflected  pulse;  directly 
on  the  opposite  side  of  the  magnetosphere  from  the  shock 
the  corresponding  L  values  were  6.44  and  7.83.  However, 
electrons  with  the  optimum  energy  and  azimuth  at  L  =8.5 
were  convected  to  L  =  2.26.  An  electron  starting  with  an 
energy  of  2.53  MeV  at  L~ 8.5  reached  a  final  L  value  of 
2.27  and  a  final  energy  of  21  MeV.  An  electron  with  an 
initial  energy  of  8.99  MeV  and  an  initial  position  of  L  =  9 
acquired  the  most  energy.  Its  final  energy  was  49.3  MeV  at 
i= 2.99.  Electrons  starting  at  lower  L  values  with  the  op¬ 
timum  initial  local  time  and  energy  will  always  penetrate 
to  lower  L  values  than  electrons  starting  at  larger  L  values 
with  the  same  magnetic  moment;  however,  no  electrons 
in  our  simulations  penetrate  below  L  —  2.  We  conclude 
that  energetic  electrons  drifting  in  phase  with  the  wave- 
field  perturbation  can  move  to  much  lower  L  values  than 
low  energy  electrons  which  track  the  magnetosphere  dis¬ 
tortion.  Figure  3  shows  the  time  evolution  of  a  ring  of 
electrons  initially  at  L  =  7  with  an  energy  of  4  MeV.  This 
figure  illustrates  the  distortion  of  the  initial  energetic  elec¬ 
tron  distribution  produced  by  the  shock.  Rings  at  other  L 


Fig.  2.  Simulated  electron  energy  spectra  at  L=  2.25,  2.45, 
2.65  and  2.85  ±A L  =  0.1  in  each  case  and  averaged  over  a 
drift  period.  Vertical  axis  is  logarithmic  in  counts  vs.  en¬ 
ergy  in  MeV,  and  spectra  are  binned  in  1  MeV  increments 
and  plotted  at  the  upper  limit  of  each  bin.  Input  simu¬ 
lation  spectrum  at  L  =  8.5  is  plotted  as  the  dot-dashed 
power  law  on  the  left. 


values  with  the  same  magnetic  moment  are  nested  within 
or  outside  the  deformed  ring  in  a  similar  manner. 

The  electron  energy  gain  is  very  sensitive  to  the  lo¬ 
cal  time.  Those  electrons  which  stay  in  phase  with  the 
reflected  pulse  longer  than  with  the  incoming  pulse,  or 
encounter  the  incoming  pulse  on  the  night  side  where  it 
is  relatively  weak,  may  lose  energy.  Figure  4  shows  the 
trajectories  of  three  electrons  starting  at  L  =8.5.  Two  of 
the  electrons  reach  a  final  L  value  of  2.5  while  the  third  es¬ 
capes  the  simulation  region  (L  >  12).  Electron  A  reaches 
an  energy  of  37.3  MeV  starting  from  an  initial  energy  of 
5.52  MeV  while  electron  B  reaches  a  final  energy  of  13.1 
MeV  starting  with  an  energy  of  1.71  MeV.  Electron  C 
starting  with  an  energy  of  5.52  MeV  escapes  the  magne¬ 
tosphere  by  interacting  more  strongly  "ith  the  reflected 
pulse  than  with  the  incoming  pulse.  Note  „hat  during  part 
of  its  trajectory  electron  A  moves  radially  inward.  This  is 
because  the  magnetic  field  of  the  inward  propagating  pulse 
cancels  the  normal  inward  gradient  of  the  dipole  magnetic 
field  and  thus  for  a  while  there  is  little  gradient  drift.  Dur¬ 
ing  this  phase  the  electron  gains  its  energy  only  from  the 
second  term  in  (2).  Electrons  like  B ,  however,  contribute 
much  more  to  the  flux  at  L  =2.5  since  the  assumed  ini¬ 
tial  energy  distribution  (IT-8)  implies  that  there  are  many 
more  such  lower  energy  electrons. 


Fig.  3.  Time  evolution  of  an  electron  ring  in  L  —  <j>  space, 
at  50  s  intervals  from  the  beginning  of  the  simulation.  The 
electrons  were  placed  initially  at  L  =  7  with  an  energy  of 
4  MeV  and  longitudinal  spacing  of  0.5°. 


Electrons  injected  at  a  given  point  at  time  t0  will  ar¬ 
rive  at  the  detector  at  time  t  =  t0  +  6/vp  where  6  is  the 
distance  along  the  drift  path  and  vp  is  the  drift  velocity 
which,  for  relativistic  electrons,  is  closely  proportional  to 
their  energy.  This  implies  that  1  jvp  =  (1  /6)(t  —  f0).  Thus 
the  points  on  a  plot  of  1/vp  versus  t  should  be  aligned 
along  a  straight  line  with  a  slope  of  1/6  and  an  intercept 
of  to  if  all  the  electrons  were  injected  at  the  same  location 
and  time.  Figure  5  is  such  a  plot  for  the  first  three  drift 
echoes  for  all  the  electrons  that  passed  0300  MLT  between 
2.46  <  L  <  2.54.  This  figure  can  also  be  interpreted  as  an 
energy-time  plot  of  the  electrons  with  the  energy  given  on 
the  right-hand  side.  Since  all  electrons  are  plotted  as  dots 
there  is  no  weighting  involved.  To  a  good  approximation, 
the  electrons  do  in  fact  fall  on  a  line  from  which  it  is  possi¬ 
ble  to  determine,  following  the  reasoning  above,  that  they 
were  injected  at  the  approximate  time  to  =  83  s  and  at  a 
position  corresponding  to  1930  MLT.  This,  however,  does 
not  correspond  to  where  the  SSC'  impacted  the  magneto¬ 
sphere  in  the  model  and  as  inferred  from  data.  Instead,  the 
simulated  perturbation  was  centered  at  1500  MLT  and  the 
electrons  were  accelerated  over  finite  local  times  and  radial 
distances.  The  straight  dashed  line  in  Figure  5  corresponds 
to  the  injection  of  all  electrons  at  1500  MLT  and  L  =  2.5. 
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Fig.  4.  Simulated  trajectories  of  three  equatorially  mir¬ 
roring  electrons  starting  at  L= 8.5.  Electrons  A ,  B,  and  C 
start  at  0900,  1500,  and  2100  MLT,  respectively. 
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consequence  we  confirm  that  the  shock  impacted  the  mag¬ 
netosphere  in  the  afternoon  sector,  that  it  compressed  the 
magnetosphere  inside  geosynchronous  orbit,  probably  to 
about  4  Re,  and  that  the  electrons  were  energized  to  about 
50  MeV  in  less  than  100  s.  In  addition  we  show  that,  al¬ 
though  the  injection  time  and  location  is  well-defined  by 
the  energy  dispersion  of  the  simulated  electrons,  this  lo¬ 
cation  does  not  correspond  to  the  actual  region  where  the 
electron  energization  took  place.  Future  work  will  include 
the  effects  of  a  variable  pulse  speed  due  to  the  interaction 
with  the  plasmasphere,  off-equatorial  electrons  and  the  ac¬ 
celeration  of  ions,  including  the  double  peak  effect  noted 
here. 
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Fig.  5.  Plot  of  l/t>o  vs.  t  for  all  particles  in  the  simulation 
which  intersect  L  =  2.5  ±0.04  at  0300  MLT.  The  three  sets 
of  points  correspond  to  the  first  three  drift  echoes.  The 
dashed  line  is  the  loci  of  points  expected  if  all  the  particles 
were  injected  simultaneously  at  1500  MLT  at  L— 2.5. 

In  contrast  to  the  local  injection  scenario,  the  lower  energy 
electrons  are  accelerated  closer  to  the  spacecraft  where  the 
wave  field  is  weaker  while  higher  energy  electrons  are  ac¬ 
celerated  farther  from  the  spacecraft  where  the  wave  field 
is  stronger.  Injections  commonly  seen  at  geosynchronous 
orbit  often  have  energy  dispersion  from  which  injection  dis¬ 
tances  can  be  calculated  [Belian  et  al.,1978].  Based  on  our 
results,  one  should  be  very  cautious  in  inferring  the  mag¬ 
netic  local  time  of  injection  or  energization. 

Another  feature  of  Figures  3  and  5  is  the  hole  in  the 
electron  flux  at  energies  between  12  MeV  and  30  MeV. 
Had  these  electrons  been  measured  by  a  detector  with  nar¬ 
row  energy  resolution,  this  gap  would  have  appeared  as 
a  double  peak  in  the  electron  flux  vs.  time  at  these  en¬ 
ergies.  Double  peaks  were  observed  in  the  ion  flux  dur¬ 
ing  this  event  [Blake  et  al.,  1992]  and  are  commonly  seen 
in  ion  drift-echo  events  at  geosynchronous  orbit  [Belian  et 
al.,1984].  The  hole  maps  to  the  space  between  A  and  B 
in  Figure  3.  Due  to  the  lack  of  electrons  at  L  >  9  in  the 
model,  empty  phase  space  extends  inside  L  =  2.5. 

Conclusions 

The  most  important  conclusion  of  our  simulation  is  that 
the  electron  injection  and  drift-echo  event  of  March  24, 
1991  can  be  understood  as  a  consequence  of  a  simple  model 
of  the  electric  and  magnetic  field  perturbations  in  the  mag¬ 
netosphere  initiated  by  the  interaction  of  an  interplanetary 
shock  with  the  magnetosphere.  While  our  model  of  the 
electric  and  magnetic  fields  is  simplified,  and  other  models 
may  produce  similar  results,  it  incorporates  enough  realis¬ 
tic  features  to  model  this  electron  drift-echo  event.  As  a 
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